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1.  Introduction. 

The  combination  of  laser  spectroscopy  and  supersonic  molecular  jet 
expansions  has  made  possible  the  study  of  a  wide  array  of  weakly  bound  van  der 
Waals  (vdW)  molecules  in  the  gas  phase.  These  clusters,  formed  in  the  jet 
expansion,  are  stable  in  the  post-expansion  region  and  can  be  studied  as 
isolated  molecules.  They  are  interesting  both  theoretically  and  experimentally 
because  of  their  unique  characteristics  such  as  low  binding  energies,  large 
intermolecular  equilibrium  distances,  and  low  frequency  intermolecular  vibra¬ 
tional  vdW  modes.  Furthermore,  the  vdW  clusters  only  slightly  perturb  the 
individual  properties  of  their  molecular  constituents.*  These  characteristics 
set  the  vdV  cluster  apart  as  a  distinct  phase  of  matter  to  be  explored  and  under 
stood . 

The  electronic-vibrational  spectroscopy  of  aromatic  molecules  like 
benzene^”4  and  s-tetrazine'*  clustered  with  various  solvents  reveals  inter¬ 
esting  information  regarding  unique  cluster  characteristics.  Specifically,  the 
studies  show  detailed  information  pertaining  to  the  intermolecular  energetics 
and  dynamics  of  cluster  interactions,  especially  in  the  area  of  the  clusters' 
low-frequency  vdV  vibrational  modes.  These  modes  are  of  considerable  interest 
since  they  represent  the  precursors  of  a  variety  of  condensed  phase  eigenstates 
such  as  phonons  in  liquids  and  solids.  They  are  also  key  factors  in  cluster 
dynamics  as  they  play  a  major  role  in  the  energy  transfer  processes  of  intra¬ 
molecular  vibrational  redistribution  (IVR)  and  vibrational  predissociation  (VP). 

Presently,  little  information  is  available  concerning  the  detailed 
energetics  of  vdV  nodes  either  experimentally  or  theoretically.  The  theoreti¬ 
cal  work  on  cluster  energetics  and  dynamics,  to  date,  treats  only  the  vdV 
16  7 

stretching  mode.  ’  *  In  these  studies,  the  clusters  are  treated  using  a 
"dumbbell"  approximation  in  which  the  intermolecular  motion  is  restricted  to  a 
stretching  mode  form.  Within  this  approximation,  the  neglect  of  bending  and/or 


torsional  vdV  inodes  seriously  limits  the  application  of  existing  theoretical 
treatments  to  IVR  and  VP  phenomena  occurring  in  molecule-molecule  type  clusters. 
The  stretching  mode  restriction  dictates  that  ovly  certain  energy  transfer  pro¬ 
cesses  can  occur,  specifically,  those  which  involve  the  vdW  Ptretch.  Experi¬ 
mental  evidence  of  the  theory's  limitation4  is  borne  out  in  the  observation  of 
energy  transfer  from  prepared  states  tangential  to  the  vdW  stretching  motion. 
Part  of  the  difficulty  of  incorporating  the  vdW  bending  and/or  torsional  inodes 
into  the  theory  is  that  no  model  has  been  developed  to  describe  either  their 
energetics  or  mode  nature. 

O 

In  previous  publications,  we  have  used  theoretical  calculations  to 
elucidate  structure  and  binding  energy  in  the  vdV  cluster  ground  state. 
Combining  these  calculations  with  experimental  observables  such  as  binding 
energies,  ionization  energies,  spectral  shifts,  relative  feature  intensities, 
and  the  appearance  of  cluster  constituent,  symmetry  forbidden  transitions  has 
aided  considerably  in  spectral  assignment  and  understanding.  The  calculations 
always  yield  results  consistent  with  experimental  observations  in  regard  to 
the  number  of  cluster  configurations  observed,  their  respective  binding 
energies,  and  their  qualitative  geometries. 

In  the  majority  of  the  cluster  spectra  studied,  vdV  vibronic  features 
are  observed  and  sometimes  assigned  based  upon  overtone  and  combination  band 
analysis.  In  other  cases,  however,  these  vibronic  features  are  numerous  and 
conplex,  making  elucidation  of  the  mode  fundamentals  difficult  and  sometimes 
impossible.  The  vdW  vibronic  feature  assignments  have  been  made  based  upon  tne 
assumptions  tnat  the  vdW  stretch  occurs  at  higher  frequency  and  with  greater 
intensity  than  vdV  bends  and  torsions.  Assigning  spectra  using  these  assump¬ 
tions  is  a  difficult  task  without  a  priori  knowledge  of  vibronic  mode  nature 
since  one  must  consider  that,  in  nonlinear  polyatomic  molecule-molecule 
clusters,  six  vdV  modes  exist  of  which  only  one  ie  a  stretching  mode.  (Ir  a 
non-linear  atom-polyatomic  molecule  cluster,  three  vdW  modes  exist). 


The  studies  reported  in  this  publication  involve  the  calculation  of  a 

complete  set  of  ground  state  vdW  vibrational  modes  for  benzene(Ar)1 , 

s-tetrazine(Ar)^ ,  benzene(CH^)1 ,  benzeneCHgO)^ ,  and  benzeneCNH^)^ . 

The  calculations  are  performed  using  a  self-consistent  pairwise  atom-atom 

a 

intermolecular  potential  developed  by  Scheraga,  et  al.  containing  general 
nonbonding  (6-12),  general  hydrogen  bonding  (10-12),  and  monopole  charge 
parameters.  The  calculated  ground  state  vdV  modes  are  compared  with  cluster 

2  7  5 

vibronic  spectra  previously  studied  in  this  and  other  laboratories. 
Consequently,  a  number  of  spectral  reassignments  are  suggested.  Vibronic 
selection  rules  governing  the  vdW  cluster  S^SQ  transition  are  also  derived 
based  upon  calculated  ground  state  cluster  mode  symmetry,  cluster  geometry  and 
experimental  observation.  Within  this  framework,  ramifications  of  Herzberg- 
Teller  vibronic  coupling  are  discussed  as  they  pertain  to  experimental 
observations  of  Franck-Condon  forbidden  transitions. 

The  calculated  vdV  vibrations  are  presented  for  all  systems  considered 
as  eigenvector  normal  modes  and  eigenvalue  energies  determined  via  normal 
coordinate  analysis  of  the  entire  vdW  cluster.^®  In  performing  the  calcula¬ 
tions,  the  high  frequency  intramolecular  vibrations  of  the  cluster  constituents 
are  assumed  to  be  completely  uncoupled  from  the  low  frequency  vdW  modes. 

Simpler  models  are  also  considered  in  studying  the  vdW  vibrations  of 
benzene(Ar)^ .  This  system  is  studied  using  four  methods:  (1)  a  Taylor  series 
expansion  of  the  intermolecular  vdV  potential  along  the  three  Cartesian  axes  in 
which  the  term  coefficients  are  related  to  the  vibrational  frequencies  in  these 
directions;  (2)  a  Morse  potential  fit  to  the  intermolecular  vdW  potential  along 
each  of  the  three  Cartesian  axes;  (5)  a  Morse  potential  fit  to  the  intercole- 

cular  vdW  potential  using  $  ■  —  ;  ’  and  (4)  a  semi-classical  energy 

Ro 

level  fit  to  the  intermolecular  vdW  potential  using  the  JWKB  method .  These 
studies  are  presented  to  show  the  consistency  of  the  cslculations  and  to  reveal 


the  advantages  and  pitfalls  of  the  models.  The  four  models  involve  treating 
the  vdW  clusters  as  simple  two  particle  systems  whose  motions  are  restricted  to 
the  principle  axes  of  the  Cartesian  coordinate  system.  The  intermolecular  vdW 
potential  surfaces  are  thereby  reduced  to  one-dimensional  potential  functions 
which  can  be  analyzed  in  a  diatomic  molecule  approximation. 

The  above  outlined  approaches  are  all  based  more  or  less  on  a  rigid 


molecule,  bound  potential  energy  well  approximation.  This  approach  would  seem 

quite  reasonable  for  the  stretching(s)  and  bending  (b  and  b  )  degrees  of 

x  y 

freedom  which  are  essentially  translations  of  the  components  of  the  vdW 

molecules  with  respect  to  one  another.  Torsional  modes  (t^,  t  ,  tz),  on 

the  other  hand,  could  in  principle  be  modeled  by  a  free/hindered  rotor 

formalism.  *  A  one  dimensional  "free  rotor"  description  has  been  applied 

to  the  symmetry  axis  (z)  torsional  motion  of  benzene  and  toluene  (CH^, 

(CD^)^  and  (CF^)1  and  compared  to  the  experimental  observations. 

Preliminary  results  suggest  that  this  approach  does  not  faithfully  reproduce 

the  experimentally  observed  spectra  for  this  "isotopic"  series  in  terms  of 

line  shapes,  intensities,  major  features,  and  the  number  of  observed 

transitions.  A  three  dimensional  "free  rotor"  model  has  also  been  applied  to 

this  problem  in  order  to  treat  all  torsional  modes  (t  ,  t  ,  t  ) 

x  y  z 

simultaneously.  Similar  difficulties  are  experienced  in  fitting  the 
experimentally  observed  spectra.  An  account  of  these  studies  will  be 
submitted  for  publication  in  the  near  future. 

The  driving  motivation  in  these  studies  is  to  answer  the  following 


questions:  (l)  if  parametric  calculations  involving  cluster  geometry  and 
binding  energy  are  consir  tent  with  experiment,  can  the  same  data  set  be 
utilized  to  calculate  intermolecular  vibrational  modes;  and  (2)  whet  are  the 


advantages  and  pitfalls  of  the  various  models  in  regard  to  the  complexity  of 
calculation,  the  approximations  made,  and  the  nature  of  the  results  obtained? 

II.  Experimental  Procedures. 

Experimental  data  pertaining  to  the  benzene(Ar)^  vdW  vibronic 
spectrum  are  obtained  employing  the  experimental  apparatus  and  procedures 
similar  to  those  used  previously  to  study  vdW  clusters.  The  benzene(Ar)1 
SQ  spectrum  is  recorded  using  a  pulsed  supersonic  molecular  jet  expansion 
in  combination  with  1 -color  time-of-f light  mass  spectroscopy  (TOFMS).  A  single 
Md  vYAC  pumped  LDS  698  dye  laser  whose  output  is  frequency  doubled  and  then 
mixed  with  the  Nd  ^/YAG  1 .064  m  fundamental  is  used  to  probe  the  6^  region 
of  the  benzene(Ar)^  cluster.  A  5%  Ar  in  He  mixture  is  placed  inline  with 
liquid  benzene  in  a  trap  maintained  at  room  temperature.  This  three  component 
gas  mixture  is  then  expanded  using  a  pulsed  nozzle  maintained  at  100  psig 
backing  pressure.  Apparatus  chamber  pressure  is  maintained  at  or  below  4  x  10 
torr  during  the  experiment. 

III.  Theoretical  Considerations. 

The  normal  coordinate  analyses  of  the  vdW  clusters  are  conducted 
employing  the  GF  methods  of  Wilson.^  These  methods  involve  solving  the 
characteristic  equation  of  5N-6  coupled  harmonic  oscillators  for  its  3N-6 


non-sero  eigenvalues  and  eigenvectors.  The  approach  is  to  treat  the  vdW 
cluster  as  a  "giant  molecule"  and  treat  both  the  intramolecular  vibrational 
modes  and  the  internolecular  vdW  modes  simultaneously.  The  intennolecular  vdW 
potential  field  used  in  the  analyses  is  expressed  in  an  internolecular 
symmetry  coordinate  system.  In  this  coordinate  system,  the  internolecular 
force  field  is  diagonal  in  the  3N  dimensional  space. 

The  cluster  constituent  intramolecular  vibrational  frequencies  are 
considerably  higher  than  those  of  the  vdW  modes.  A  reasonable  approximation 
in  this  context  then  is  to  assume  that  the  intramolecular  modes  are  completely 
uncoupled  from  the  low  frequency  vdW  modes.  Thus,  the  intramolecular  modes 
are  -aken  to  be  those  of  the  cluster  constituents.  The  constituent  force 
fields  are  generated  using  the  centrel  force  approximation11  including  out- 
of-plane  motion  terms.  Since  the  vdW  modes  arise  from  the  restriction  of 
cluster  constituent  translations  and  rotations  and  the  central  force  approxi¬ 
mation  adequately  reproduces  these  degrees  of  freedom,  the  necessary  uncoupling 
of  the  intramolecular  modes  and  the  intennolecular  vdV  modes  is  maintained  along 
with  providing  adequate  vdV  mode  calculational  results.  Other,  more  sophisti¬ 
cated  force  field  approximations  are  tested  in  the  calculations;  the  central 
force  approximation  is  determined  to  be  adequate  for  calculations  of  the  vdV 
modes.  The  only  restriction  which  applies  is  that  the  intramolecular  field 
yields  mode  eigenvalues  in  the  proper  frequency  regions. 

Within  the  central  force  approximation,  the  force  field  contains  only 
diagonal  terms  in  the  internuclear  symmetry  coordinate  system.  This  diagonal 
force  field  provides  a  simple  and  convenient  mesne  of  using  approximate  force 
constants  for  intramolecular  motion  in  the  calculations.  The  intramolecular 
force  constants  chosen  are  those  pertaining  to  general  functional  group 
stretches  and  bends.11 

The  most  convenient  choice  for  a  coordinate  system  as  the  working 
bans  for  matrix  diagonalisstion  is  the  Cartesian  system.  This  coordinate 


system  is  chosen  since  it  is  the  system  used  in  the  present  cluster  configur- 
stion  calculations,  and  more  importantly,  it  is  the  coordinate  system  in  which 
the  G  matrix  is  diagonal  and  obvious. 


i 


In  order  to  combine  the  intenr.olecular  and  intramolecular  force  fields 
algebraically,  the  force  fields  are  transformed  into  the  Cartesian  coordinate 
system.  The  transformation  yields  two  3N  dimensional  P  matrices.  The  intra¬ 
molecular  matrix  is  designated  as  Fq  and  the  intermolecular  matrix  is  designated 

as  F'.  The  F  matrix  consists  of  two  diagonal  blocks  containing  the  coordinates 
—  — o 

of  cluster  constituent  intramolecular  motion.  The  F‘  matrix  contains  two  off- 
diagonal  blocks  and  diagonal  entries  corresponding  to  the  "perturbations"  yield¬ 
ing  the  vdV  motion.  Adding  these  two  matrices  results  in  the  "giant  molecule" 

F  matrix  of  order  3N.  This  matrix  is  left-multiplied  by  the  G-1  matrix  and 
numerically  diagonalized  with  the  eigenvalues  and  eigenvectors  being  determined 
in  the  usual  fashion.  Upon  diagonalization,  the  eigenvalues  and  eigenvectors  of 
intramolecular  motion  are  identified  along  with  those  corresponding  to  cluster 
translation  and  rotation.  These  latter  zero  energy  modes  are  discarded  and  the 
remaining  nodes  are  the  eigenvalues  and  eigenvectors  of  the  ground  state  vdW 
modes. 


The  intermolecular  force  constants  used  in  the  normal  coordinate 
analysis  are  generated  from  the  intercolecular  vdV  potential  by  making  a 


harmonic  approximation.  Within  this  approximation,  the  force  constant  is 

9 


simply  the  second  derivative  of  the  potential  function,- 


l< »  » 


n  t 

* 

|2L 

k: 

_  _  i 

i  t  . 

►  :  j-: 

I'  •«“ 

*  1 

i  i  • 

1*21. 

.  lIL 

!«* 

I',- 

• « 
t 

;  • 

JJJ.C  s,*, 


Or 


tj  kl 


l''V-a 


t  i  r  » 

;  •  nr 


The  potential  function  contains  a  general  nonbonded  potential  (NB),  a 
monopole  electrostatic  potential  (ME),  and  a  general  hydrogen  bonded  potential 
(HB)  in  a  Lennard-Jones  (6-12-1-10-12)  form.  The  total  intermolecular  inter¬ 
action  is  taken  as  a  sum  of  pairwise  aton-atom  interactions  over  all  the  atoms 
of  each  cluster  constituent.  The  r. .'s  are  the  atoD-stom  distances  between 

atom  i  on  constituent  k  and  atom  j  on  constituent  1.  The  r.  ,'s  represent,  the 

—  J 

symmetry  coordinates  in  which  the  intermolecular  force  field  is  diagonal.  The 

second  derivative  of  eq.  (l)  with  respect  to  r. .  gives  the  force  constant  of 

^  J 

the  atom-atom  interaction  as. 
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in  which  the  K.  .'a  are  elements  of  the  F'  matrix  expressed  in  intermolecular 
symmetry  coordinates.  These  terms  are  evaluated  at  the  equilibrium  configura¬ 
tion  of  the  cluster  assuming  that  the  cluster  constituents  are  frozen  with 
regard  to  intramolecular  motion. 

Three  additional  models  are  employed  to  study  the  benzene(Ar)1  vdV 
cluster  as  a  test  case.  In  these  models,  the  vdW  cluster  is  assumed  to  be  a 
"diatomic  molecule"  in  the  sense  that  the  system  is  considered  to  be  composed 
of  two  particles,  the  benzene  molecule  (solute)  and  the  argon  atom  (solvent). 
The  vdW  modes  are  assured  to  arise  from  restricted  1 -dimensional  motion  of 
the  benzene  molecule  relative  to  the  argon  atom.  Only  the  atom-molecule 
benzene(Ar)1  cluster  will  be  considered  using  these  models  since  its  vdV 


modes  are  easily  characterized  by  simple  translations  along  any  of  tne  three 
Cartesian  axes.  The  vdW  stretching  mode  is  taken  as  motion  restricted  to  the 


1 -dimensional  translation  moving  the  cluster  constituents  apart  in  opposite 
directions.  The  vdV  bending  inodes  are  considered  to  be  notions  restricted  to 
1 -dimensional  translations  which  move  the  cluster  constituents  parallel  to  one 
another  in  opposite  directions. 

In  the  first  model  considered,  1 -dimensional  potential  curves  are 
mapped  out  by  translating  the  solvent  molecule  relative  to  the  solute  atom  in 
one  of  tne  three  Cartesian  directions.  The  intermolecular  vdW  potential  is 
assumed  to  be  represented  by  a  Taylor  series  expansion  about  the  equilibrium 
intermolecular  distance,  Bq  in  the  form, 


KR)  •  l(R  i 
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The  expansion  coefficients  are  evaluated  by  a  polynomial  fit  to  the  potential 
curves  taking  the  displacement  vector  R  as  the  independent  variable  and  U(R) 
as  the  dependent  variable.  The  second-order  polynomial  fit  coefficient 
determines  the  effective  harmonic  force  constant  governing  the  frequency  of 

i  ? 

bound  state  motion.  The  energy  of  thi6  motion  is  given  by 
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The  third  end  fourth  order  polynomial  fit  coefficients  represent  an  enharmonic 

correction  term  to  the  energy.  From  perturbation  theory,  the  correction  can 

1 2 

be  written  to  first  order  as 
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Higher  order  terms  are  neglected  in  the  enharmonic  corrections  since  u  » 

e 

ueXe  >>  other  corrections.  A  polynomial  least-squares  fit  to  tenth  order  in 

R  is  determined  to  be  sufficient  to  reproduce  faithfully  the  1 -dimensional 

potential  curves  generated  via  translation.  In  passing  we  note  that  the  Taylor 

series  expansion  could  also  be  evaluated  directly  by  taking  successive 

derivatives  of  eq.  (l).  This  nay  be  the  method  of  choice  if  one  is  only 

interested  in  the  lowest  order  terms  in  evaluating  eq.  (4)  or  (5>),  or  if  one 

is  interested  in  a  more  "exact"  reproduction  of  the  potential  curves  in  the 

vicinity  of  the  dissociation  limit. 

The  second  model  considered  involves  fitting  the  1 -dimensional 

potential  energy  curves,  derived  in  the  same  manner  as  previously  discussed, 

1 S 

to  a  Morse  function  of  the  form, 


L(R)  -  D  <e”2r(R"Ro>  .  2e~MR"Rc'V 


The  entrgv  level*  and  enharmonic  corrections  are  evaluated  using 


The  last  model  used  involves  an  energy  level  fit  based  upon  the  semi- 
classical  JWKE  method.^  In  this  model,  the  energy  levels  are  determined 
from  the  quantization  of  the  action  integral  according  to  the  Bohr-Sommerfeld 
restrictions.  The  governing  equation  is 

I  -  (2d)”**  /  [£  -  ueff(R)^  dR  *  h(v  +  'i)*  (9) 

Energy  level  determination  is  accomplished  by  numerical  integration  of  (9) 
taking  ueff(R)  as  the  1 -dimensional  intermolecular  vdW  potential  in  a  tenth 
order  polynomial  form.  The  path  of  integration  is  taken  over  one  complete 
notion  cycle  with  the  boundary  conditions  established  by  the  intermolecular 
potential  at  a  specific  energy  E. 

IV.  Results. 

A.  Benzene(Ar)j . 

Figure  1  and  Table  I  present  the  benzene(Ar)1  vdW  cluster 

spectrum  recorded  in  the  region  between  38561  cr.”1  and  30710  cn-1  using 

1 -color  TOFKS.  The  cluster  is  "red  shifted”  by  21  cm  1  with  respect 

to  the  benzene  b\  The  bathvchromic  shift  is  indicative  of  the  greet er 

binding  energy  in  the  cluster  S,  state  relative  to  the  cluster  S  state. 

i  o 

Three  vibronic  features  are  observed  to  the  blue  of  the  cluster  61 .  No 

0 

features  are  observed  in  the  symmetry  "forbidden”  benzene  0°  region; 
therefore,  the  cluster  must  have  at  least  a  three-fold  axis  of  symmetry. 


Figure  2  and  Table  II  contain  the  celculational  results  of  the  ground 

state  configuration  and  vdW  modes  of  benzene(Ar)^ .  Configurational 

calculations  yield  a  single  geometry  of  minimum  energy  for  the  cluster 

possessing  C^v  symmetry,  figure  2.  In  this  geometry,  the  argon  atom  lies 

3.44  X  above  the  benzene  molecular  plane  along  the  z  (6-fold)  axis.  The 

ground  state  cluster  binding  energy  is  calculated  at  287  cm  1  which  makes 

the  excited  state  binding  energy  308  cm"1 .  The  calculated  intermolecular 

distance  of  3*44  X  compares  well  with  that  of  3*4  X  obtained  from  the 
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rotational  analysis,  adding  independent  proof  to  the  adequacy  of  the 
calculations  in  predicting  detailed  information  on  cluster  structure. 

The  eigenvalue  energies  from  the  normal  coordinate  analysis  of 
benzene(Ar)1  are  40  cm-1  and  11  cm”1  for  the  vdW  stretch  and  vdW  bends, 
respectively.  The  eigenvector  normal  modes,  figure  2,  reveal  that  the  vdW 
stretch  entails  purely  perpendicular  motion  of  the  argon  atom  relative  to  the 
benzene  molecular  plane.  Furthermore,  the  calculations  reveal  that  the  2-fold 
degenerate  vdV  bending  mode  involves  some  combination  of  motion  parallel  to 
the  benzene  molecular  plane.  Both  of  these  eigenvector  results  are  consistent 
with  group  theoretical  arguments  as  expected. 

The  "diatomic  molecule"  model  calculations  yield  three  sets  of  vdV 
mode  energies.  The  average  mode  frequencies  are  40  cm"1  for  the  vdW  stretch, 
«2(a,)  and  10  cm”1  for  the  vdW  bends  b^(e^).  Figures  3  and  4  show  the 
details  of  the  1 -dimensional  potential  curve  mappings  and  the  results  of  the 
model  calculations,  Note  that  the  z-directior.  potential  curve  modeling  the 
vdW  stretch  looks  surprisingly  similar  in  form  to  that  of  a  typical  diatomic 
molecule.  All  models  yield  both  adequate  potential  curve  fits  and  consistent 
vibrational  energy  level  structures. 

B.  a-Tetrasine^Ar)^ . 

Figure  3  presents  the  results  of  the  ground  state  configuration  and 
vdW  mode  analysis  of  s-tetrazine( Ar)^ .  Only  the  normal  coordinate  analysis 


vibrational  calculation  ie  presented  since  this  method  yields  the  most  inform¬ 
ative  results  for  our  purposes  and  the  consistency  between  the  eigenvalue/ 
eigenvector  results  and  the  "diatomic  molecule”  results  has  already  been  shown 
for  the  benzene(Ar)^  case.  Configurational  calculations  yield  a  single  cluster 

geometry  of  minimum  energy  possessing  C^v  symmetry,  figure  5*  In  this  geo- 

o 

metry,  the  argon  atom  lies  3.45  A  above  the  s-tetrazine  molecular  plane  along 

the  z  (2-fold)  axis.  The  calculated  ground  state  cluster  binding  energy  is 

295  cm”1 .  The  calculated  intermolecular  distance  of  3«45  £  compares  extremely 

o 

well  with  the  intermolecular  distance  of  3*45  A  obtained  from  the  rotational 

analysis.  The  calculated  binding  energy  of  295  cm”1  also  lies  within  the 

"  -15 

experimental  limits  of  254  <  D  <  332  cm  .  Again,  the  calculations  and 

o 

experiment  are  in  exact  agreement. 

The  normal  coordinate  analysis  eigenvalues  are  41 .0  cm”1  for  the  vdW 

stretching  mode  s  (a,),  9  cm-1  for  the  vdW  bending  mode  b  (b,),  and  12 
zi  y  i 

cm”1  for  the  vdW  bending  mode  The  eigenvector  normal  modes, 

figure  5,  show  that  the  vdW  stretch  is  restricted  to  motion  perpendicular  to 
the  s-tetrazine  molecular  plane  while  the  vdW  bends  are  restricted  to  motion 
parallel  to  the  molecular  plane.  As  in  the  benzeneCAr)^  analysis,  these  results 
are  consistent  with  group  theoretical  arguments. 

C.  Benzene(CH  ), . 

The  ground  state  configuration  and  vdW  eigenvalues/eigenvectors 

for  benzene(CH4)1  are  shown  in  figure  6.  The  results  presented  for  the 

geometry  and  binding  energy  of  benzene(CK  )  are  in  good  agreement  with 
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previous  reports  from  this  laboratory  using  an  exponential-six  and  Lernard- 

2 

Jones  potential  form.  In  this  geometry,  figure  6,  the  methane  center-of- 
o 

mass  lies  at  3*47  A  above  the  benzene  molecular  plane  on  the  principle  s  (3-fold) 
axis.  The  cluster  ground  state  binding  energy  is  540  cm”1 .  The  norr.al  coor¬ 
dinate  analysis  reveals  six  vdV  vibrations,  two  being  2- fold  degenerate.  The 


ground  state  vibrational  energies  are  82  cr."1  for  the  vdV  stretch  s  (a.), 

2  1 

16  cm”1  for  the  vdV  bends  b  (e),  and  28  cm”1  and  89  cm”1  for  the  vdW 

xy 

torsions  t  (a0)  and  t  (e),  respectively.  The  eigenvector  normal  modes, 
s  t  xy 

figure  6,  transform  as  the  translational  and  rotational  representations  of  the 
Cjv  point  group,  as  indicated.  The  vdW  stretching  mode  transforms  as  the 
translation  of  the  cluster  constituents  away  from  one  another  along  the  z  (3- 
fold)  axis.  The  vdV  bending  modes  transform  as  some  combination  of  cluster 
constituent  translations  in  opposite  directions  perpendicular  to  the  3-fold 
axis  in  the  xy  plane.  One  vdW  torsion  mode  transforms  as  s  rotation  of  the 
cluster  constituents  about  the  z  (3-fold)  axis  in  opposite  directions.  The 
remaining  two  vdW  torsions  transform  as  rotations  about  orthogonal  axes  per¬ 
pendicular  to  the  3-fold  axis. 

D.  Benzene(H20)1 . 

The  calculated  benzene(H20)1  geometry  used  in  the  normal 

3 

coordinate  analysis  is  similar  to  that  calculated  previously.  Only  one 

minimum  energy  configuration,  which  has  a  binding  energy  of  504  cm”1,  is 

found.  As  shown  in  figure  7,  the  cluster  geometry  possesses  symmetry 

o 

with  the  H^O  center-of-mass  located  3-15  A  above  the  benzene  molecular  plane. 

Six  ground  state  vdW  vibrations  sre  calculated  for  the  C  cluster 

8 

geometry.  Their  corresponding  eigenvalues  and  eigenvectors  are  shown  in 
figure  7.  The  six  vdW  modes  consist  of  a  vdW  stretch  at  159  cm”1,  two  vdW 
bends  st  14  cm”1  and  18  cm"1,  and  three  vdW  torsions  at  40  cm”1,  50  cm”1, 
and  15b  cm’1.  The  eigenvector  normal  mode3  transform  as  the  translational  and 
rotational  representations  of  the  C  point  group:  the  vdV  stretch  transforms 

B 

as  a  z  translation;  the  vdVi  bends  transform  as  x  and  y  translations;  and  the 

vdV  torsions  transform  as  R  ,  R  ,  and  R  rotations. 

x  y  i 

E.  Benzene(NH^)^ . 

Configurational  calculations  on  the  benzenetKK^  cluster 
reveal  two  cinimur  energy  geom®triee  similar  to  those  obtained  previously.^ 


One  cluster  geometry  possesses  C^v  symmetry  with  a  binding  energy  of  711 
cm”1  while  the  other  possesses  C  symmetry  and  s  binding  energy  of  608 

B 

cm-1 .  In  the  C^v  cluster,  figure  8,  the  NH^  center-of-msss  is  located 
o 

3.23  A  above  the  benzene  molecular  plane  along  the  z  (3-fold)  axis.  In  the 

0 

Cs  cluster,  figure  9,  the  NH^  center-of-mass  is  located  3*29  A  above  the 
benzene  molecular  plane. 

Using  the  potential  surfaces  generated  from  these  two  configurations, 

six  ground  state  vdW  vibrations  are  calculated  for  each  geometry.  Their 

corresponding  eigenvalues  and  eigenvectors  are  shown  in  figures  8  and  9*  The 

C3v  cluster  ground  state  normal  modes  transform  in  the  same  manner  as  those 

of  benzene ( CH .  In  the  C  cluster,  the  ground  state  normal  modes  trans- 
4  1  8 

form  similar  to  those  of  benzeneCl^O^  . 

V.  Discussion. 

To  compare  the  calculated  and  experimentally  observed  vdW  modes,  we 

assume  that  the  intermolecular  potential  surface  of  the  cluster  is  identical 

in  both  Sq  and  S1  electronic  states.  This  assumption  is  justifiable  if  one 

considers  that  cluster  fluorescence  excitation  and  dispersed  emission  spectra 

4 

are  similar  for  the  vdW  vibronic  transitions.  Furthermore,  the  small 

spectral  shifts  of  the  chromophore  and  the  week  intensity  of  the  vdW  modes 

signify  only  small  changes  in  cluster  binding  energy.  This  is  probably 

indicative  of  only  slight  variations  of  the  potential  surface  between  S  and 

0 

S1  electronic  states. 

Comparisons  between  calculation  and  experiment  are  made  using  group 
theoretical  arguments  based  on  the  selection  rules  governing  the  vibronic 
transitions.  Specifically,  transition  moment  matrix  elements  are  qualitatively 
analyzed  using  the  crude  adiabatic  approximation  for  which  the  vibrational  code 
dependence  on  the  electronic  wave  function  is  explicit.  In  this  case,  a  stand¬ 
ard  Herzberg-Teller  (H?)  expansion  and  adiabatic  wave  functiors  are  used14 


with  the  electronic  wave  function  vibrational  rode  dependence  truncated  at 
second  order. 

Using  this  expansion,  two  unique  types  of  spectra  can  be  generated. 
First,  one  could  consider  that  the  vdV  nodes  do  not  participate  in  the 
vibronic  coupling  schene  (Case  I)  and  that  they  merely  enter  into  the 
expansion  as  an  additional  scalar  product  (overlap  integral).  This  argument 
dictates  that  only  totally  symmetric  Franck-Condon  progressions  and  combination 
bands  are  spectroscopically  observed.  Furthermore,  the  intensities  of  these 
features  are  solely  derived  from  the  cluster  chromophore  vibronic  mode  with 
which  they  are  in  combination.  Alternately,  one  could  consider  the  vdW  modes 
to  be  capable  of  vibronic  coupling  (Case  II),  In  this  case,  they  enter  into 
the  HT  transition  moment  equation  in  the  aSme  manner  as  other  vibronically 
active  modes.  The  operator  responsible  for  these  transitions  would  be  of  the 
form 

/  »2u  — \  e  Q6qvdw 

V  3Q*3qvdv  /  Q6  -  Qe.  qvdw  *  qvdw 

in  which  is  the  cluster  vibrational  mode  and  qydW  is  *  specific 
cluster  vdV  node.  This  argument  allows  the  possibility  of  observing 
noctotally  symmetric  vdV  fundamentals  with  "borrowed"  intensity  due  to 
interelectronic  state  mixing. 

In  the  individual  cluster  discussions,  the  above  two  cases  are 
considered  in  order  to  assign  ar.d  understand  the  observed  cluster  vibrunic 
spectra. 

A.  Bensene(Ar)1 . 

Comparison  of  the  calculated  ground  state  vdW  vibrations  of 
bencer.e(Ar)1  and  the  experimental  vibronic  spectrum,  figure  1,  Table  I, 


and  Table  II,  shows  that  vibronic  assignments  cen  be  made  based  upon  calcu¬ 
lations  for  all  models  analyzed.  Considering  the  transition  moment  matrix 
elements  for  the  benzene(Ar)1  S1-*-So  transition  under  C^v  symmetry,  one 
should  expect  to  observe  totally  symmetric  combination  bands  of  the  vdW 

stretch  built  upon  the  cluster  61.  The  selection  rule  for  these  combination 

0 

bands  is  dv  ■  0,  H,  ♦ 2 ,  ...  Also,  one  should  expect  to  observe  nontotally 
symmetric  vdW  bend  combinations  with  the  6^  with  the  selection  rule  being 
dv  *  0,  +2,  . . •  These  selection  rules  hold  for  both  Case  I  and  Case  II  type 

spectra  and  imply  that  the  vdW  modes  do  not  enter  into  the  vibronic  intensity 
borrowing  mechanism. 

The  calculated  vdW  stretching  mode  at  40  cm-1  compares  quite  well 
with  the  experimental  vibronic  feature  at  39*7  cm-1  to  the  blue  of  the 
benzene(Ar)j  6^.  Thus,  this  feature  is  assigned  to  the  benzene(Ar)1  vdW 
stretch/cluster  61  combination  band  61  sfa,)1  based  upon  the  dv  ■  0, 

♦1,  ^2  ...  selection  rule. 

For  the  vdW  bending  modes,  only  odd  overtones  are  expected  to  be 

observed  as  pointed  out  above.  The  experimentally  observed  features  at  30.9 

cm  *  and  61.8  cm  1  to  the  blue  of  the  cluster  6^  correspond  to  overtone 

features  of  the  vdW  bends;  using  the  dv  ■  0,  * 2 ,  *4  ...  selection  rule,  the 

former  feature  is  the  first  overtone  of  the  vdW  bends  and  the  latter  is  the 

third.  Considering  the  30.9  cm-1  feature  as  the  first  overtone  places  the 

symmetry  forbidden  bend  fundamental  at  about  15*5  cm”1.  This  energy  lies 

close  to  the  calculated  2-fold  degenerate  ground  state  vdV  bend  at  11  cm”1. 

Thus,  these  two  spectral  features  are  assigned  to  the  benzene(Ar)j  vdW  bends 

1  12 

first  and  third  overtone/cluster  6  combination  bands,  6  b  (e,)  and 

o  o  xy  i  o 
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Ramifications  of  Herzberg-Teller  vibronic  coupling  in  the  benzene(Ar) 
cluster  are  obvious.  Fron  the  derived  vibronic  selection  rules  and  experi- 
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Rental  observation,  the  benzene(Ar)^  cluster  spectrum  is  best  assigned  based 

upon  benzene  Herzberg-Teller  coupling  (i.e.,  the  6^  feature  is  allowed)  with 

vdW  totally  symnetric  modes  and  combinations  forming  short,  weak  Franck-Condon 

progressions  built  upon  the  intense  benzene  transition.  In  addition,  the  cal- 

1 5 

culated  and  observed  cluster  symmetry  is  verified  by  the  vdW  vibronic 

structure. 

B.  a-Tetrazine(Ar)j . 

The  selection  rules  governing  the  s-tetrazine(Ar)  vdW  vibronic 

transitions  under  C?v  symmetry  arise  from  Case  I  Franck-Condon  arguments. 

The  totally  symmetric  vdW  stretch  should  be  observed  to  the  blue  of  the 

cluster  0°  following  a  Av  -  0,  +1 ,  >2,  ...  selection  rule.  The  vdW  bends 
o  —  — 

should  only  be  observed  in  odd  overtones  (Av  *  0,  ♦ 2 ,  M  ...)  built  on  the 
allowed  s-tetrazine  0®  transition.  As  in  the  benzene(Ar)1  cluster,  no  Case 
II  distinction  can  be  made  for  s-tetrazine(Ar)1  and  thus  no  vdW  Herzberg- 
Teller  vibronic  coupling  is  expected.  Unfortunately,  a  complete  experimental 
spectrum  showing  the  details  of  the  s-tetrazine(Ar)^  vdW  modes  is  not,  as 
yet,  available.  The  only  information  in  this  regard  is  the  identification  of 

the  vdW  stretching  mode  at  44  cr”1  to  the  blue  of  the  cluster  0®  by  Levy 

5  -1-1 

et  al.  Other  vdW  features  at  66  cm  and  106  cm  to  the  blue  of  the 

cluster  0®  are  observed,  but  they  are  neither  assigned  nor  are  their 

spectra  published. 

Figure  10  and  Table  III  compare  the  calculated  ground  state  vdW  mode 
frequencies  with  those  observed  in  the  s-tetrazine( Ar),  1 B,  ♦  *A  vdW 
vibronic  apectrur.  The  experimentally  assigned  vdW  stretch  et  44  cm-1  to 
the  blue  of  the  cluster  0®  corresponds  to  the  41  cm”1  calculated  ground 
state  stretch.  This  motion,  like  that  in  the  benzene(Ar)j  cluster,  involves 
perpendicular  motion  of  the  argon  ator.  relative  to  the  s-tetrazine  molecular 


17 


plane.  Based  upon  ground  state  calculations,  the  feature  66  cm'1  to  the  blue 

of  the  cluster  0°  probably  corresponds  to  a  vdW  stretch/vdW  bend  overtone 

12  12 

combination  band,  s  (a4)  b  (b0)  or  s(a.)  b  (b, )  ;  the  feature 

£102^0  i  o  y  »  o 

at  108  cm”^  probably  corresponds  to  the  next  allowed  bend  overtone/stretch 

combination  band  s(a, b  (b0)l  or  s(a,)^  b  (b. )^.  No  speculation 

■  o  x  £  o  i  o  y  i  o 

with  regard  to  which  bend  is  responsible  for  the  features  observed  will  be  made 
since  no  spectra  of  sufficient  sensitivity  are  available  for  analysis. 

C.  Benzene (CH^)^ 

In  a  previous  analysis  of  the  benzene(CH^)^  cluster,  three 
major  vdW  vibronic  features  were  reported  in  the  cluster  6^  region:  these 
features  were  assigned  to  a  bend  fundamental  (27.3  cm'1),  a  stretch  funda¬ 
mental  (32.3  cm-1),  and  a  stretch  overtone  (51*4  cm-^).  The  assignments 
were  made  based  upon  the  assumptions  described  in  the  Introduction. 

Considering  the  H7  transition  moment  matrix  elements  and  assuming  a 
Case  I  type  spectrum,  the  selection  rules  for  the  benzene(CH . vdV  vibronic 
transitions  are  av  *  0,  +1 ,  *2,  ...  for  the  vdW  stretch  and  Av-  0,  ♦ 2 ,  +4,  ... 

for  the  vdW  bends  and  torsions  when  in  combination  with  the  cluster  61 . 

0 

The  selection  rules  involved  in  Case  II  can  be  viewed  in  two  ways.  The 

selection  rules  can  be  derived  using  the  calculated  cluster  symmetry  of  C^y, 

or  they  can  be  derived  by  considering  the  cluster  symmetry  as  C^.  The  latter 

situation  arises  since  the  metnane  center-of-mass  is  calculated  at  3*47  A  above 

the  benzene  molecular  plane.  At  this  distance,  the  methane  could  be  viewed  as 

a  sphere  above  the  benzene  molecular  plane  and,  hence,  the  use  of  the  C^y 

point  group  to  represent  the  vdW  mode  symmetries  could  be  warranted. 

In  Cj^  symmetry  Case  II,  the  selection  rules  for  vdW  mode  combinations 

with  the  cluster  61  are  av  *  0*  +1,  ♦2,  ...  for  all  six  vdW  modes.  If  C, 

0  —  —  ?  v 

is  the  correct  cluster  physical  symmetry  and  H?  vdV  coupling  exists,  all  modes 
can  be  observed  in  the  cluster  6^  region.  In  C^y  symmetry.  Case  II,  the  selec¬ 


ts 


-  -A  Vi  j v. 


tion  rules  for  the  vdV  mode  combinations  with  the  cluster  6^  are  4v  ■  0,  ♦! ,  <2 

...  for  the  vdW  stretch  a  (a^)  and  torsion  tz(a2)  and  Av  «  0,  <2,  *t ,  ... 

for  the  2-fold  degenerate  vdW  bends  b^Ce^)  and  torsions  t  (e^).  In  this 

approximate  high  symmetry,  only  the  cluster  61  and  vdW  mode  t  are  capable  of 

z 

vibronic  coupling  for  the  6^  transition. 

Comparison  of  the  experimental  6^  vibronic  spectrum  of  benzene(CH^  )1 

and  the  calculated  ground  state  vdV  vibrations  is  shown  in  figure  12  and  Table 

IV.  The  observed  feature  at  27.3  cm”1  to  the  blue  of  the  cluster  61  corre- 

o 

sponds  to  the  t  torsion  calculated  at  28  cm*1.  Thus,  this  feature  and  its 
z 

observed  overtones  at  51*4  cm-1  and  73*5  cm-1  are  reassigned  to  61  vdW  torsion 

0 
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combination  bands  6  t  (a„)  ,  6  t  (a~)  ,  and  6  t(a0)  ,  using 

O  Z  d  0  O  Z  c.  0  0  Z  c  0 

the  Cjv  Case  II  or  C^v  Case  II  Av  *  0,  ,  *2,  ...  selection  rule.  The  Av 

-  0,  *1 ,  *2,  ...  selection  rule  suggests  that  for  the  benzeneCCH^  system 
vdW  mode  Herzberg-Teller  vibronic  coupling  is  an  important  component  of  the  over¬ 
all  intensity  mechanism.  In  the  present  case,  the  occurrence  of  the  nontotally 

symmetric  t  torsion  progression  implies  that  the  vdW  modes  are  vibronically 
z 

active  and  that  assuming  them  to  be  nonparticipants  in  the  coupling  mechanism 

oversimplifies  the  physics  necessary  to  explain  the  cluster’s  spectroscopy. 

The  intense  feature  at  32*3  cm*1  to  the  blue  of  the  cluster  61 

o 

corresponds  to  the  first  overtone  of  the  2-fold  degenerate  vdV  bending  mode, 
12  -1 

6  b  (e)  ,  calculated  at  16  cm  •  This  identification  is  based  upon  the 
o  xy  o 

observation  of  a  feature  at  64.6  cm”1  which  corresponds  to  the  third  overtone 
of  the  bends.  Furthermore,  weak  intensity  features  are  observed  at  about  16.1 
cm*1  and  48.4  cn"1  which  could  <'orrespond  to  the  vdW  bend  fundamentals  and 
second  overtones.  Tne  observation  of  tnese  features  adds  proof  to  the  arguments 
suggesting  that  the  vdV  modes  are,  at  least,  minor  participants  in  the  Herzberg- 
Teller  vibronic  coupling  scheme.  The  cluster  symmetry  and,  hence,  spectroscopy 
are  thereby  also  best  described  using  the  calculated  C^y  point  group  in  con- 


junction  with  Herzberg-Teller  coupling  rather  then  the  approximate  Cfev  point 

group.  Based  upon  this,  the  features  at  32*3  cm“^  and  64.6  cm”^  are  reassigned 

11  2  1 

to  vdW  bend  overtone  combinations  with  the  cluster  6  ;  6  b  (e)  and  6 

o  o  xy  o  o 

b^y(e)^  using  the  selection  rule  Av  *  0,  ,  * 2 ,  ....  The  features 


are  assigned  in  Table  IV  and  figure  11.  The  weak  features  at  16.1  cm  and 

48.4  cm-1  are  assigned  to  61  b  (e)1  and  61  b„,  (e)^.  Based  upon  the 

o  xy  o  o  xy  o 

relative  intensities  displayed  in  the  spectrum,  vdW  vibronic  coupling  is  an 
important  factor  in  the  intensity  of  this  progression  involving  the  e  symmetry 
bending  modes. 

Neither  the  vdW  stretch  ®z(«^)  nor  the  2-fold  degenerate  torsions 

t^y(e)  are  observed  in  the  6^  spectrum.  This  fcould  be  due  to  poor  Franck- 

Condon  factors  for  these  vibronic  transitions  since  they  are  both  calculated  to 

be  at  relatively  high  energies  (ca.  82  cm-1  and  89  cm“\  respectively). 

Moreover,  these  modes  could  be  participating  in  VP  since  the  total  energy  6^ 

a  (a.)1  or  61 t  (e)1  is  close  to  that  of,  if  not  above,  the  cluster's 
z  1  o  o  xy  o 

S1  binding  energy.  Both  the  vdW  stretch  and  t^y  torsions  involve  motion  per¬ 
pendicular  to  the  benzene  molecular  plane:  this  motion  could  couple  well  to  the 
VP  process. 


D.  Benzene(H^O)1 . 

The  Sj  ♦  SQ  vibronic  spectrun  of  benzene(HpO)1  previously 

observed  in  this  laboratory^  possesses  two  unique  spectral  regions  located 

around  the  cluster  0°  and  6^  containing  vdW  vibronic  features.  No  vdW 

oo 

vibronic  assignments  were  made  in  either  region  and  no  correlation  between  the 
regions  was  suggested. 

Examination  of  the  H?  transition  moment  matrix  elements  using  C 

8 

symmetry  and  Case  I  considerations  leads  to  the  selection  rules  iv  -  0,  H, 

22,  ...  for  the  vdW  s  stretch,  b  bend,  and  t  torsion  and  Av  «  0,  ♦ 2 , 

t  y  i  — 

♦4,  ...  for  the  b  bend,  t  and  t  torsions.  Under  Case  II  arguments, 

“  *  7  1 


the  selection  rule  is  jv  •  u,  ♦!,  *2,  ... 
modes  are  capable  of  vibronic  coupling. 


for  all  aix  vdV  modes:  all 


Vibronic  spectra  of  benzene(H-O),  in  both  the  cluster  0°  and  61 

ci  0  0 

regions  are  reproduced  in  figure  12  along  with  the  calculated  ground  state  vdW 

mode  energies.  The  observed  feature  at  5.2  cm”1  (4.8  cm”1)  to  the  blue  of 

the  cluster  0°  (61)  transition  corresponds  to  the  vdW  b  bend  fundamental 
oo  x 

calculated  at  14  cm”1.  The  observed  feature  at  16.2  cm”1  (15-8  cm”1)  to 

the  blue  of  the  cluster  0°  (61 )  transition  corresponds  to  the  calculated 

o  o 

totally  symmetric  vdW  b  bend  at  18  cm”1.  Additionally,  the  observed  features 

y 

at  34.6  cm”1  (34.6  cm”1)  and  49  cm”1  (48.4  cm”1)  to  the  blue  of  the 
cluster  0°  (b1 )  transition  are  associated  with  the  t(a")  and  t  (a')  torsion 

0  0  Z  X 

fundamentals  calculated  at  40  cm”1  and  50  cm”1 ,  respectively.  The  occurrence 
of  the  nontotally  symmetric  fundamentals  implies  that  the  a v  *  0,  +_2,  *4,  ... 
selection  rule  for  Case  I  in  which  the  vdW  modes  are  not  vibronically  coupled 
is  violated.  The  violation  suggests  that  the  vdW  modes  participate  in  the 
vibronic  coupling  scheme  and  that  the  4v  *  0,  ^1 ,  ♦ 2 ,  ...  selection  rule 
should  apply  to  all  six  vdW  modes  (Case  II).  Based  upon  this,  the  spectra  are 
best  assigned,  Table  V  and  figure  12,  using  both  nontotally  and  totally 
symmetric  vdW  progressions. 

Assigning  the  benzenetH^O)^  spectra  using  the  CB  point  group  repre¬ 
sentations  corroborates  the  calculated  cluster  geometry.  Treating  the  cluster 
in  approximate  high  symmetries,  such  as  C^v,  leads  to  selection  rules  which 
are  clearly  violated  when  applied  to  spectral  observation  and  assignment.  Spe¬ 
cifically,  the  bend  and  torsion  fundamentals  are  forbidden  under  these  higher 
symmetry  approximations.  Furthermore,  the  spectral  assignments  using  C 

8 

symmetry  arguments  suggest  that  the  water  constituent  is  likely  located  above 
tne  benzene  molecular  plane. 

Neither  the  vdW  a  stretch  nor  the  t  torsion  are  observed  in  either 

*  y 

spectral  region.  This,  as  in  the  benzenefCH^)^  case,  probably  results 
from  poor  Franck-Condon  factors  for  these  particular  modes;  they  are  both 


calculated  to  be  at  relatively  high  energies  (ca.  159  cm-1  and  156  cn"  , 
respectively).  Moreover,  in  the  6^  region,  these  modes  could  be  participating 
in  VP  since  the  totsl  energy  of  the  system  at  these  levels  is  close  to,  if  not 
above,  the  cluster  binding  energy  (c.a.  •  500  cm-1).  The  decrease  in  the 
hypsochromic  shift  and  the  shift  of  intensity  maximum  in  the  vdW  manifold  in 
going  from  the  cluster  0®  to  the  6^  may  also  be  indicative  of  a  VP  process. 

E.  Benzene(NH^)^ . 

Benzene(NH^)^  clusters  recently  observed  in  thi9  laboratory^ 

yield  spectra  in  both  the  cluster  0®  and  6^  regions.  Two  cluster 

geometries  are  calculated  for  the  system,  one  possessing  C  symmetry  and  the 

other  possessing  symmetry.  Prom  symmetry  arguments  the  C0  symmetry 

cluster  is  the  only  contributor  to  the  0°  spectrum  while  both  cluster 

geometries  contribute  to  the  6^  spectrum.  Neither  of  these  spectra  were 

analyzed  nor  assigned  in  the  initial  observation  since  they  are  bo  complicated. 

They  were  merely  presented  as  an  indication  of  the  notion  that  cluster  vibronic 

spectra  can  sometimes  be  very  extensive  and  congested. 

The  benzene(NKj)^  Cg  symmetry  cluster  follows  the  same  vdV  vibronic 

selection  rules  as  derived  for  the  benzene(H20)1  cluster.  In  this  symmetry, 

the  selection  rules  for  Case  I  are  Av  «  0,  ,  ^2,  ...  for  the  vdV  stretch, 

b  bend,  and  t  torsion  and  Av  ■  0,  j^2,  ^4,  ...  for  the  b  bend,  t  and 
y  x  x  y 

t  torsions.  In  Case  II  the  selection  rule  is  Av  ■  0,  ♦1 ,  *2,  ...  for  all 
z  —  — 

six  vdV  modes.  Furthermore,  the  calculated  cluster  geometries  are  qualitatively 
eimilar.  Hence,  their  0°  spectra  should  be  qualitatively  similar:  this  is 
borne  out  in  both  experimental  results  and  vibrational  mode  calculations. 

The  benzene(NHj)j  C^y  symmetry  cluster  follows  the  same  vibronic 
selection  rules  as  presented  for  ber.zene(CH, ), .  Here  the  selection  rules 
are  either  av  -  0,  ♦  1 ,  *2,  ...  for  the  vdV  stretch  and  Av  ■  0,  <2, 

^4,  ...  for  the  vdV  bends  and  torsions  (Case  I)  or  av  ■  0,  *1,  *2,  ... 

for  all  an  vdV  modes  (Case  II). 


r 


The  calculated  ground  atate  vibrational  energies  snd  the  obeerved 


cluster  0°  and  fe^  vibronic  spectra  are  compared  in  figure  1?.  Due  to  the 


complex  nature  of  the  spectra  which  possibly  results  from  hindered 


rotation  within  the  cluster  and  hot  bands,  only  tentative  assignments  of  the 
vdW  mode  progressions  are  made.  The  tentative  assignments  are  based  upon  both 
vibrational  mode  calculations  and  upon  inference  from  the  benzene(H^O) 1  and 
benzene(CH^)1  cluster  spectra. 


In  the  0  spectrum,  the  most  intense  low  energy  feature  is  assigned 


to  the  origin  of  the  cluster's  Sq  transition.  The  smaller  intensity 

features  to  the  red  of  this  feature  are  thus  hot  bands  which  yield  the 


sequence  structure  in  the  0°  region.  The  observed  feature  at  8.8  cm-1  to 


the  blue  of  the  cluster  Cr  corresponds  to  the  nontotally  symmetric  vdV  b^ 


-1 


bend  fundamental  calculated  at  15  cm  .  Additionally,  the  observed  features 


at  15.0  cm-1,  45*1  cm”1,  54.2  cm"1,  99-6  cm”1,  and  127.7  cm”1 


correspond  to  the  calculated  b  bend  (21  cm’1),  t  torsion  (44  cm”1), 

y  z 


t^  torsion  (48  cm  '),  stretch  (112  cm  1),  and  t  torsion  (125  cm”1). 


The  occurrence  of  the  nontotally  symmetric  fundamentals  suggests  that  the  vdW 
modes  participate  in  vibronic  coupling  (Case  II)  and  that  the  Av  *  0,  *1, 

<2,  ...  selection  rule  should  apply  to  all  six  vdV  modes.  Using  this 
selection  rule,  the  spectrum  is  best  assigned,  Table  VI,  using  both  nontotally 


and  totally  symmetric  vdV  progressions. 


1 


In  the  spectrum,  the  most  intense  low  energy  feature  is  assigned 


to  the  cluster  6^  vibronic  origin.  This  assignment  bears  a  resemblance  to 


the  benzeneO^O)^  6^  spectrum  in  the  respect  that  the  6^  features  for  both 

clusters  are  red  shifted  relative  to  that  observed  et  the  S,  S  origins  (ca. 

1  o 

25  cm’1  for  benzenefNH^J^  and  35  cm’1  for  benzenefH^Oj^ ). 


23 


Table  VI  presents  the  C  symmetry  geometry  assignments  in  the  61 

8  0 

region.  The  assignments  are  made  using  the  vdV  fundamentals  identified  in  the 

o£  spectrum  and  &v  *  *1 ,  +2,  ...  selection  rule  for  all  six  vdW  modes  (Case  II). 

The  contribution  of  the  symmetry  benzene(NHj) ^  cluster  to  the 

6^  spectrum  is  observed  starting  with  the  feature  at  72.9  cm  1  to  the  blue 

of  the  C  cluster  6* .  Assignment  of  this  feature  to  the  C,  cluster  61 
so  5v  o 

is  based  upon  the  observation  that  no  intense  features  are  seen  in  either  the 

benzene(NH^)1  0°  or  the  benzenelH^O)^  6^  spectra  at  this  energy.  The 

observed  feature  at  15  cm”1  to  the  blue  of  the  C^v  cluster  6^  corresponds 

to  the  2-fold  degenerate  vdW  b  bends  at  19  cm”1.  The  observed  feature  at 

53  cm”1  to  the  blue  of  the  C,.  cluster  61  is  associated  with  the  calcu- 

5v  o 

lated  t  torsion  at  44  cm”1.  As  in  the  benzene(CH^)^  case,  the  observation 

of  the  b  and  t  fundamentals  implies  that  the  vdW  modes  participate  in  the 
xy  z 

vibronic  coupling  mechanism  (Case  II)  and  that  the  £v  *  ,  ±2,  ...  vibronic 

selection  rule  applies  for  all  vdW  modes.  Based  upon  this,  the  best  assignments 

for  the  C^v  cluster  vdV  vibronics  are  tabulated  in  Table  VI. 

F.  Morse  Model. 

Ro 

Calculations  of  the  ground  state  vdW  modes  of  benzene(Ar)1  are 

16  7 

also  conducted  using  a  model  described  by  Jortner  et  al.  '  '  This  model 
contains  a  main  feature  which,  at  first  glance,  appears  to  make  it  generally 
applicable  in  predicting  vdV  stretching  mode  frequencies  from  calculated 
cluster  binding  energies  and  equilibrium  intermolecular  distances.  However, 
the  model  turns  out  to  be  incorrect  in  this  application.  The  model  employs  a 
simple  relationship  between  the  cluster  equilibrium  intermolecular  distance 
H0  and  the  Morse  potential  parameter  f  .  The  relation  is  derived  by  equating 
the  second  derivative  of  eq.  (t>)  with  the  second  derivative  of  a  single  term 
(b-12)  potential  function.  This  derivation  yields  the  relation  6  "  ^ 

o 

which  can  be  used  to  predict  the  vdV  stretching  mode  frequency  by  diatomic 


Korse  fit  methods. 


Substituting  8*  |  and  the  calculated  ground  state  cluster  binding 
o 

energy  into  eq.  (7)  yields  a  vdVf  stretching  mode  energy  of  47  cn  for 

benzene(Ar)^  employing  Scheraga's  potential  parameters.  Even  though  the 

calculated  vibrational  energy  is  qualitatively  correct,  figures  14  and  3  show 

that  the  8  *  g  relation  results  in  an  inadequate  reproduction  of  the 
o 

potential  curve,  especially  in  the  critical  region  near  the  equilibrium  inter- 
molecular  distance.  Additional  calculations  using  the  potential  data  described 
in  ref.  1  and  the  above  approximate  model  also  result  in  an  inadequate  potential 
curve  reproduction,  figure  15«  In  this  case,  the  vdV  stretching  mode  is  cal¬ 
culated  at  55  cm”^. 

The  failure  of  the  model  under  both  data  sets  suggests  that  the  model 

itself  is  inadequate  in  this  application.  The  model  fails  in  this  application 

since  the  6*  |  relation  holds  exactly  only  in  Bystems  in  which  the  molecule- 
Ho 

atom  Lennard-Jones  parameters  have  been  determined  directly.  Applying  this 
approximate  model  to  a  case  in  which  the  potential  function  is  represented  by 
pairwise  atom-atom  potentials  does  not  take  into  account  the  differing  contri¬ 
butions  of  each  interaction  to  the  potential  energy  and  the  equilibrium  inter- 
molecular  distance.  In  this  more  complicated  situation,  no  analytic 
relationship  between  8  and  exists. 

Comparing  the  vdV  stretching  mode  energies  calculated  using  Jortner'a 
and  Scheraga’s  data  sets  suggests  that  a  significant  difference  exists  between 
the  two  data  sets.  Benzene(Ar)^  configurational  calculations  using  Jortner’s 

data  aet  yield  a  single  cluster  geometry  of  C^v  symmetry  with  the  argon  atom 
o 

located  3*5  A  above  the  benzene  molecular  plane.  In  this  case,  the  ground  atete 
cluster  binding  energy  is  395  cm”1.  This  binding  energy  is  10b  cr.”1  greater 
than  that  calculated  using  Scheraga's  data  set.  The  binding,  energy  of  ?b7 
cc."1  calculated  using  Scheraga's  data  aet  is  probably  more  accurate  since  the 
benzeneCAr)^  ground  state  binding  energy  should  be  very  similar  to  that  cal- 


culated  for  s-tetrazine(Ar)^ .  in  the  latter  case,  the  calculated  ground  atate 

binding  energy  of  295  cfi"1  compares  well  with  that  observed  experimentally, 

11  -15 

254  <  D  <  332  cm  .  Furthermore,  generating  Lennard- Jones  parameters 
o 

using  the  data  set  of  ref.  1  yields  a  binding  energy  of  359  cm”1  for 
s-tetrazine(Ar)^  which  is  clearly  not  as  accurate  as  the  binding  energy 
reported  in  this  work.  The  difficulty  is  due  to  the  consolidation  of  atom-atom 
parameters  from  different  data  sets.  In  this  respect,  great  caution  must  be 
taken  when  consolidating  parameters  since  each  parameter  set  is,  in  general. 


only  self-consistent  and  may  have  no  meaning  when  combined  with  parameters  from 
other  sets. 


VI.  Conclusions. 


TOFKS  studies  have  been  employed  to  determine  the  general  geometry  and 
symmetry  of  vdV  clusters  in  the  gas  phase.  Through  computer  modeling,  a 
correlation  between  the  details  of  the  cluster  geometry  and  Bpectral  features 


has  been  demonstrated.  Specifically,  the  parametric  calculations  yield  useful 
information  regarding  cluster  geometry,  binding  energy,  and  the  vdW 
vibrations.  These  calculated  results  are  consistent  with  experiment  and  serve 
as  predictive  and  analytic  tool  which  can  be  used  to  elucidate  and  understand 
the  details  of  vdV  cluster  energetics. 

Of  the  several  models  considered  in  studying  the  intemolecular  vdV 


modes,  simple  diatomic  approximations  yield  adequate  results  when  applied  to 
atom-molecule  clusters.  On  the  other  hand,  for  molecule-molecule  clusters 


normal  coordinate  analysie  is  essential.  Tne  normal  coordinate  analysis  is 
especially  useful  for  analyzing  systems  which  have  little  or  no  symmetry  since 
no  a  priori  knowledge  of  vdW  mode  nature  ia  necessary  to  generate  potential 


energy  aurface  mappings. 

Reessigneenta  of  and  assignments  to  cluster  vdV  modes  have  been  cade 
based  upon  the  knowledge  gained  froc  calculation,  from  comparison  of 


TO 


■"1  Ir  « 


calculation  and  experiment,  several  conclusions  result.  First,  the  actual 
excited  state  normal  mode  vdW  frequencies  are  well  fit  by  the  calculated 
ground  state  cluster  potential.  This  conclusion,  though  not  surprising,  gives 
independent  proof  of  the  invariance  of  the  weak  vdW  potential  between  ground 
and  excited  electronic  states  as  well  as  providing  a  means  of  using  ground 
state  vdW  vibrational  structure  to  predict  vdW  vibronic  structure.  Second,  in 
the  majority  of  the  clusters  analyzed,  the  observed  vdW  vibrations  are  those 
involving  bending  and  torsional  motions  parallel  to  the  aromatic  ir-system. 
Furthermore,  these  modes  are,  in  general,  quite  active  in  the  Herzberg-Teller 
vibronic  coupling  mechanism  and  significant  interelectronic  state  mixing 
obtains.  Third,  vdW  motions  which  penetrate  the  aromatic  w-system  have  high 
frequencies  and  are  only  observed  in  systems  in  which  the  Franck-Condon 
factors  and  binding  energies  are  favorable.  Finally,  the  observed  vibronic 
structure  supports  the  calculated  cluster  geometry  in  all  cases. 
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TABLE  I 


vdW  spectral  features  in  toenzeneUr^  61  region  and  calculated  ground  state 
vdW  modes  (refer  to  fig.  1).  0 


Energy  Relative  to  Calculated  Ground  a 

Cluster  6^  (cm-1)  State  Energy  (cm-1) 

0  (38587.6) 

11  <v 

30.9 

39.7  40  ( s  ) 

z 

61.8 


Assignment 


o 


61*,  2 

o  xyo 


,1  1 
6  8 
0  zo 


6!b  1 
o  xyO 


a)  vdV  mode  representations  as  per  fig.  2. 


TABLE  II 


Calculated  ground  atate  vdW  node  energies  for  benzenefAr)^ . 


Model 

s  (cm”1 )  8 
z 

b  (cm”1  ) 

xy 

Taylor  Series 

40.88  (1.51) 

9.80  (.03) 

Morse  Fit 

39.47  (1.3b) 

10.54  (.10) 

JVKB 

40.05  (1.43) 

9-71  (.02) 

Normal  Coordinate 
Analysis 

40.0 

11.0 

a)  Values  in  parentheses  are  first  order  anharmonicity  corrections  calculated 
from  "diatomic  molecule"  models.  vdW  mode  representations  as  per  fig.  2. 


TABLE  III 


vdV  spectral  features  in  s-tetrazine(Ar)1  0°  region  (s-tetrazine  1  B^u  ■*-  1 A  ) 
and  calculated  ground  state  vdW  inodes  (refer  to  fig.  10). 


Energy  Relative  to 

Cluster  0°  (cm"' )  C 
o 

0  (18104.9) 


44  a 
66 

108 


Calculated  Ground 

State  Energy  (cm"1)  D  Assignment 


9  (by) 

12  (\) 

41  (sz) 


o 


zo 


8  ’b  2 

ZO  zo 


8  ’b  4 
zo  zo 


or  s 


or  s 


zo  yo 

»b  4 
zo  yo 


a)  Observed  and  assigned  in  ref.  5« 

b)  vdV  node  representations  as  per  fig.  5» 

c)  From  ref.  5* 


TABLE  IV 


vdW  spectral  features  in  benzene(CH^)1  6^  region  and  calculated  ground  a 
vdV  modes  (refer  to  fig.  11). 


Energy  Relative 

to  Cluster  61  (cm”')  a 
o 

Calculated  Ground 

State  Energy  (cm”1) 

Assignment  * 

0  (38567.6) 

61 

0 

16.1 

16  (b  ) 

6\  1 

xy 

o  xyo 

27.3 

28  (t  ) 
z 

6^  t  1 
o  zo 

32.3 

6\  2 
o  xyo 

48.4 

6^b  5 

o  xyo 

51.4 

61t  2 

0  zo 

64.6 

6^b  4 

o  xyo 

73-5 

82  (s  ) 

6\  5 

0  zo 

z 

89  (t^) 

a)  From  ref.  2  and  unpublished  spectra. 

b)  vdV  mode  representations  as  per  fig.  6. 


TABLE  V 


vdW  spectral  features  in  benzene(HpO) 1  0°  and  61  regions  and  calculated 
ground  state  vdW  modes  (refer  to  fig.1 12^.  0 


Energy  Relative  to 

Cluster  0°  ^cm"') 

or  6  (cm~  ) 
o 


Calculated  Ground 

State  Energy  (cm“  )  Assignment 


0  (38168.6) 

o  o 
o 

5.2 

14  (b  ) 

b 1 

X 

xo 

16.2 

18  (b J 

b  T 

y 

yo 

21 .4 

bxib,o 

25 

b  2b 1 
xo  yo 

31.9 

b  2 
yo 

34.6 

40  (t  ) 

t  1 

z 

zo 

39.4 

45.7 

b  5 

yo 

49 

50  (t  ) 

t  1 

X 

yo 

50.8 

tzobyo 

67.4 

t  ’b  1 
yo  yo 

70.3 

t  2 
zo 

99.5 

t  2 
yo 

103.5 

t  5 
zo 

• 

156  (ty) 

- 

159  (.,) 

0  (58655*4) 


TABLE  V  (Continued) 


15.8 

18  (by) 

25.5 

30.3 

34.6 

40  (tz) 

39-6 

48.4 

50  (tx) 

60.5 

97.8 

101 .0 

156  (t  ) 

•7 

159  (.,) 

a)  From  ref.  3* 

b)  vdW  mode  representations  as  per  fig.  7. 


61  b> 


1 

o'yo 


6\  2b  1 
o  xo  yo 


6’b  2 
o  yo 


1 


61 1 
o  zo 


1.  1 


6*t  'b 
o  zo  xo 


61 1  1 
o  yo 


6'b  4 
o  yo 


61 1  2 
o  yo 


61 1  3 
o  zo 


TABLE  VI 


vdV  spectral  features  in  benzene(NH,)1  0°  and  61  regions  and  calculated 
ground  state  vdW  modes  (refer  to  fig.  13?. 


Energy  Relative  to 
Cluster  0° 
or  6  ?c«n* ' ) 


Calculated  Ground  .  . 

State  Energy  (cm  )  Assignment 


-39-2 

-36.7 

-33-4 

-31-3 

-29-6 


v  o.  ° 

i«byl 

y.  O.  1 

by2bxo 


-26.7 


-19.2 
-  6.7 


-  3-8 


0  (38021.1) 


15  (bx) 

21  U>y) 


b  f  b  f 

xl  yl 


24.2 


b  ’b  1 
xo  yo 


29-2 


b 

yo 

b  2b  1 
xo  yo 


44  (tg> 


48  (tx) 


b  ’b  2 
xo  yo 


TABLE  VI  (Continued) 


45.1 

44  (t  ) 

b’t  1 

48.7 

z 

0  zo 

53.7 

61 1  1 

0  zo  zo 

55-3 

48  (t  ) 

61 1  1 

z 

o  zo 

58.6 

61 1  S  1 

o  zo  yo 

61 .0 

61 t  ’b  2 

0  zo  zo 

64.7 

61 1  \  1 

0  zo  zo 

68.8 

72.9  [0  (38587.6)] 

‘J<V 

74.6 

«'t  ’b  2 

0  zo  ZO 

79-5 

e't  S  S  1 

o  zo  zo  yo 

83-3 

61 1  S  2 

o  zo  yo 

87.9 

19  (b  ) 

6’b  1 (C, 

o  zyo'  3v 

98.7 

112  (s  ) 

,1  1 

6  B 

z 

o  zo 

100.4 

61 1  \  1 
o  zo  zo 

103-9 

6^  2(C,  ) 

o  *yov  3vy 

107.1 

A1  K  1 

6  s_  b 
o  zo  zo 

111.9 

61 t  \  ’b  1 
o  zo  zo  zo 

114.0 

,1  K  1 

6  B  b 
o  zo  yo 

123-4 

,1  1v  1 

6  s  b  b 
o  zo  zo  yo 

125-9 

44  (tt) 

- 

97  (st) 

- 

- 

152 

- 

a)  From  ref.  3* 

b)  vdV  node  representations  as  per  fig.  8  and  9* 

c)  C-  cluster  61  contributions  tabulated  relative  to  CB  cluster  61 

vlbronic  origin  as  in  fig.  13.  0 


? I CURE  CAPTIONS 


Figure  1 


Figure  2 


Figure  3 


Figure  4 


Mass  selective  *  SQ  spectrum  and  calculated  ground  state 

vdW  modes  of  benzene(Ar)^ .  Energy  scale  is  relative  to 

benzene(Ar),  6^  transition  (3B587.6  cm*^).  Nozzle  backing 
i  o 

conditions:  P0  *  100  psig,  TQ  *  300  k.  Peak  positions  and 
assignments  as  per  Table  I  and  fig.  2. 

Calculated  ground  state  minimum  energy  configuration  (a)  and 

eigenvalue/eigenvector  vdW  modes  (b)-(d)  for  benzene(Ar)1 . 

Cluster  symmetry  is  C^v  with  an  equilibrium  intermolecular 
o 

distance  of  3*45  A.  Eigenvectors  are  normalized  and  displayed 
at  2x  magnification  (2  X  total  displacement). 

z-direction  (vdV  stretch)  potential  energy  mapping  of 
benzenetAr)^ .  Coordinate  system  is  as  shown  in  fig.  2.  Trans¬ 
lation  is  along  z  axis  with  x  and  y  coordinates  at  equilibrium 

intermolecular  distance  values.  Translation  is  displayed 

o 

relative  to  equilibrium  intermolecular  distance,  3*43  A.  (6-12) 
potential  energy  mapping  is  represented  by  o;  Taylor  series 
expansion  and  energy  levels  are  represented  by  —  ;  Morse  fit 
potential  energy  curve  is  represented  by  - — .  Vibrational  mode 
constants  as  per  Table  II. 

x(y)-direction  (2-fold  degenerate  vdV  bend)  potential  energy 
mapping  of  benzene(Ar)^ .  Coordinate  system  is  as  shown  in 
fig.  2.  Translation  is  along  x(y)  axis  with  z  and  y(x) 
coordinates  at  equilibrium  intermolecular  distance  values. 


Translation  is  displayed  relative  to  equilibrium  interaolecular 
distance  3*45  A.  6-12  potential  energy  mapping  ia  represented 
by  o;  Taylor  aeriea  expansion  and  energy  levels  are  represented 
by  — ;  Morse  fit  potential  energy  curve  is  represented  by 
— .  Vibrational  mode  constants  as  per  Table  II. 

Figure  3  Calculated  ground  state  minimum  energy  configuration  (a)  and 

eigenvalue/eigenvector  vdV  modes  (b)-(d)  for  a-tetrazine(Ar)1 . 

Cluster  symmetry  is  C2V  with  an  equilibrium  intermolecular 
o 

distance  of  3-45  A.  Eigenvectors  are  normalized  and  displayed 
at  2x  magnification  (2  A  total  displacement). 

Figure  6  Calculated  ground  state  minimum  energy  configuration  (a)  and 

eigenvalue/eigenvector  vdV  modes  (b)-(g)  for  benzene(CH . ). . 

4  • 

Cluster  symmetry  is  C^y  with  an  equilibrium  intermolecular 
o 

distance  of  3«47  A.  Eigenvectors  are  normalized  and  displayed 

o 

at  2x  magnification  (2  A  total  displacement). 

Figure  7  Calculated  ground  state  minimum  energy  configuration  (a)  and 

eigenvalue/eigenvector  vdV  modes  (b)-(g)  for  benzene(H20)1 . 

Cluster  symmetry  is  Ca  with  an  equilibrium  intermolecular 
o 

distance  of  3*13  A.  Eigenvectors  are  normalized  and  displayed 
at  2x  magnification  (2  A  total  displacement). 

Figure  8  Calculated  ground  state  minimum  energy  configuration  (a)  and 

eigenvalue/eigenvector  vdV  modes  (b)-(g)  for  benzene(HHj)^ . 

Cluster  symmetry  is  Cjy  with  an  equilibrium  intermolecular 
o 

distance  of  3*23  A.  Eigenvectors  are  normalized  and  displayed 

o 

at  2x  magnification  (2  A  total  displacement). 


Figure  9 


Calculated  ground  atate  minimum  energy  configuration  (a)  and 


eigenvalue/eigenvector  vdV  modee  (b)-(g)  for  benzene(NK^)^ . 

Cluater  symmetry  is  C  with  an  equilibrium  intermolecular 

0 

O 

distance  of  3-29  A.  Eigenvectors  are  normalized  and  displayed 
at  2x  magnification  (2  A  total  displacement). 


Figure  10 


Schematic  *B,  *A  spectrum  (from  ref.  5)  and  calculated 
3u  g 

ground  state  vdW  modes  of  s-tetrazine(Ar)^ .  Energy  scale  is 

relative  to  s-tetrazine(Ar).  0°  transition  (18104.9  cm”1). 

1  o 

Relative  feature  intensities  are  not  shown.  Feature  positions 
and  assignments  as  per  Table  III  and  fig.  5»  vdW  bends  are 


represented  by  b^  in  schematic  spectrum  (see  text  for  explana¬ 
tion). 


Figure  11  Mass  selective  ♦  SQ  spectrum  (from  ref.  2)  and  calculated 

ground  state  vdV  modes  of  benzene(CH4)1 .  Energy  scale  is 

relative  to  benzene(CH  ) .  61  transition  (38567.6  cm  1). 

4  1  o 

Feature  positions  and  assignments  as  per  Table  IV  and  fig.  6. 


Figure  12  Maas  selective  S1  SQ  spectra  (from  ref.  3)  and  calculated 

ground  state  vdW  modes  of  benzene ( H^0)1 .  Energy  scale  is 
relative  to  benzene(H_0),  0°  and  61  transitions  (38168.6  cm”1 

elO  0 

and  38655*4  cm’1,  respectively).  Feature  positions  and  assign¬ 
ments  as  per  Table  V  and  fig.  7.  a  and  t  vdV  modes  are  not 

*  y 

shown. 


Figure  13  Maas  selective  ♦  Se  spectrum  (from  ref.  3)  end  calculated 

ground  atate  vdV  modes  of  benzene(NH^)^ .  Energy  acale  is 


relative  to  benzene(NH-).  0°  and  61  transitions  for  C 

>10  0  8 

cluster  (38021.1  cm'1  and  38514.7  cm"1,  respectively). 

Feature  positions  and  assignments  as  per  Table  VI  and  fig.  8 

and  9*  C-  cluster  s  and  t  .  modes  are  not  shown. 

z  xy 

Figure  14  z-direction  (vdW  stretch)  potential  energy  mapping  of 

benzene(Ar)^  and  $  *  |  Morse  fit.  Translation  is  relative 

°  o 

to  equilibrium  intennolecular  distance,  3.45  A.  (6-12)  potential 

energy  mapping  using  data  set  from  ref.  9  is  represented  by  o; 

Taylor  series  expansion  and  energy  levels  are  represented 

toy  — ;  B  *  "  Morse  fit  potential  energy  curve  is  represented 
Ko 

toy  - .  Taylor  series  vibrational  mode  constants  as  per  Table 

II.  B  *  t  Moree  fit  vibrational  constants  are  m  *  47.22 
Ho  e 

cm”1  in  x  *  1.94  cm”1. 

©  © 

Figure  15  z-direction  (vdV  stretch)  potential  energy  mapping  of 

benzene(Ar)1  and  6  •  |  Morse  fit.  Translation  is  relative 

°  o 

to  equilibrium  intennolecular  distance,  3*5  A.  (6-12)  potential 
energy  mapping  using  data  set  from  ref.  1  is  represented  by  o; 
Taylor  aeries  expansion  and  energy  levels  are  represented  by 
—  with  ui  »  46.48  cm”1  and  o.  v  ■  1.81  cm”1 ;  B  *  4 

O  6 

Morse  fit  potential  energy  curve  is  represented  by  — .  B  ■  ^ 

R 

i  ° 

Morse  fit  vibrational  mode  constants  are  u,  ■  54.84  cm  t  > 

©  ©  © 
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